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Abstract 

Galvanic coupling of small-area (three-junction) flux qubits, using shared large Josephson junc- 
tions, has been shown to yield appreciable interaction strengths in a flexible design, which does not 
compromise the junctions' intrinsic good coherence properties. For an introduction, I recapitulate 
an elementary derivation of the coupling strength, which is subsequently generalized to the case of 
tunable coupling for a current-biased shared junction. While the ability to vary coupling constants 
by, say, 20% would be useful in experiments, sign-tunability (implying switchability) is highly 
preferable for several quantum-computing paradigms. This note sketches two ideas: a "crossbar" 
design with competing ferro- and antiferromagnetic current-biased tunable couplings, and a "me- 
diated" one involving an extra loop between the qubits. The latter is a variation on proposals 
for tunable capacitive coupling of charge qubits, and tunable inductive coupling of large-area flux 
qubits. 
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The design philosophy of the single three-junction (3JJ) qubit is to obtain a bistable 
system without relying on magnetic energy, by using multiple Josephson junctions p|. This 
revives as old idea for using multistable multi-junction loops and generalizes it to the 
quantum regime, in the presence of an arbitrary (but typically nearly half-integer) flux bias. 
As shown by Levitov et al. and apparently rediscovered by Butcher pj, the coupling 
of several such qubits can be implemented analogously, using large Josephson junctions 
inserted into the shared legs of adjacent qubits. In this way, an appreciable coupling strength 
should be easier to achieve than using inductive coupling, especially since the 3JJ loops 
typically have a small area (if anything, for a circuit of this type the design challenge will 
be to avoid the coupling strength being too large). In Levitov et a/.'s preprint, only the 
order of magnitude of the coupling strength is estimated (correctly). In Butcher's report 
(Section 3.5.1), it is stated that the Hamiltonians for Josephson and inductive coupling are 
identical. However, the former is not given let alone derived, while the correct form for 
the latter [H[ has eluded the Delft group (overestimation by a factor two) both in Butcher's 
thesis [Eq. (3.3)] and for some time thereafter 0. 

In this note, therefore, first the coupling strength will be derived, recapitulating a result 
meanwhile published in |2j. Both the result and the derivation turn out to be very similar 
to the inductive case, confirming Butcher's statement. However, the calculation is decidedly 
easier, since the leading answer is found by studying the classical potential for vanishing 
inductance (as opposed to expanding the full Hamiltonian to first order in the inductances). 




3 X '6 
FIG. 1: Two Josephson-coupled 3JJ qubits. 



The circuit is shown in Fig. ^ The relevant part of the Hamiltonian reads simply 
Uj = Yl^=o ~ Ej cos0j. For vanishing inductances, this is subject to flux quantization 
0i + 02 + 03 + 0o = 0x an d 04 + 05 + 06 — 0o = X (note the minus sign in the latter 
relation), where 4>^ b = 27r<3>"' V^o are external flux biases in phase units. We focus on the 
simplest case (fi® ,b = ir, £4,3,4,6 = E, E 2 ^ = aE (~ < a < 1), Eq ^> E, so that the potential 
becomes 

Uj((f>) = -E o cos0 o + E[— cos0! - cos 03 + acos(0i+0 3 +0o) 

— COS (f>4 — COS 06 + a COs(04+06— 0o)] • (1) 

The wells, corresponding to the classical qubit states, are partly characterized by OqUj = 
diUj = d 3 Uj = d^Uj = 8qUj = 0. Like for a single 3JJ qubit, the relevant solutions 
corresponding to actual minima are readily verified to have 

01 = 03 , 04 = 06 , (2) 
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leaving one with three nontrivial extremum equations. 
Two solutions are readily found as 



,FM =0 ™ = FM = ±arccos ( J_) . (3) 



U 3 ^) = -E -^ + 2a^E. (4) 

Clearly, these correspond to ferromagnetic (FM) configurations, in which the sense of su- 
percurrent rotation is the same in both loops. In this symmetric device, the currents in the 
central leg therefore cancel, so there is no phase difference across this leg and the remaining 
junctions are in the same state as for degenerately biased free 3JJ qubits. 
On the other hand, for the antiferromagnetic (AF) configurations, one has 



if = -4>f 



(5) 



The remaining equations 

E sin 0f = 2E sin < F , (6a) 
sin (j)f = a sin(20f +0f ) (6b) 

have to be solved perturbatively in E/E . Since the lowest order corresponds to free qubits, 
one has 

iAf i_ ______ ( l \ , A 1 ) E , <ntfr?/r?\2 



f = ±arccos(^- j + tf>— + OftE/EoY] . (7) 



Substitution into Eq. ((Baj) readily yields 



<t>f = ±2^1 - ^ + O{(E/E f] = + O[(E/E f] , (8) 

in terms of the classical equilibrium persistent current I p (H = 1). Subsequently, Eq. (}6b"j) 
leads to 0^ = ±[(2— 4a 2 )/(4a 2 — l)]I p /2eE. However, this cancels in the expansion (around 
a potential minimum) of Eq. ((TJ, and one finds 

W AF ) = W FM ) - + O[(E/E ) 2 } , (9) 

an AF coupling equivalent to a mutual inductance M e s = l/4e 2 i?o — precisely the Josephson 
inductance of the coupling junction. It should be feasible to manufacture devices with M eS 
ranging from typical magnetic values to values corresponding to a dimensionless coupling of 
order one. 

The analogy to the magnetic case [5| is complete: in Eq. Q, the energy increase 
(over the FM state and to leading order) |j5 o (0o ) 2 = I p /2e 2 E in the central junc- 
tion is overcompensated by the decrease in energy of the "qubit" junctions 2^ F ^ x [/j min , 
where d^UJ 1111 = =F/ p /2e. In fact, the coupling's AF character can be understood with- 
out doing the detailed expansion (and, hence, without requiring E to be large). Consider 
4>' = (0 FM , 0™, 0™, — FM , — FM ); in this perhaps counter-intuitive AF state, the two loop 
currents have opposite senses, yet there is no phase drop over the central junction. However, 
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this merely means that the state is non-stationary (the T-shaped islands get charged), but 
it is a well-defined point in the potential landscape, with i7j (</>') = Uj(<fi FM ) by Eq. (JTJ). 
Therefore, one necessarily has Uj(<p AF ) < Uj(4> FM ) for the AF minimum state, as is also 
seen numerically. 



For added flexibility, let us consider current-biasing the central leg |8[. This can be 
described by the potential 

U ((f)) = -E cos o - 2T0o - E[cos 4>i + cos 3 + a cos(0"-0i-03-0o) 

+ COS 4 + COS 6 + a COs(0^-04-06+0o)] , (10) 

where the phase frustrations will be specified shortly. The bias current is J x = 2eX; it may 
be large compared to the loop currents but shouldn't exceed the critical one, so \X\ < E . 
It may look asymmetric to couple the bias to the middle leg only, but the phases in the 
different legs are not independent. Pending a dynamic analysis of the full Hamiltonian, here 
we simply use (jlU|) to see what it predicts. Again minimizing w.r.t. the phases, one sees that 
(J2J) still applies. As expected, the central leg has to carry most of the bias, viz., 

0o = a«siii(|;) + ^Jr + 4 2) |j + O[(E/E ) 3 } ; (11) 

it is most consistent to retain O when evaluating the (large) first two terms of (fTUJl to 
0(E/Eq), although its actual value cancels. The need to have a stable solution unfortunately 
limits us to |0o| < n/2 in (fTTj) . 

To leading order, the outer arms should behave like degenerate 3JJ qubits, so we have to 
compensate for the contribution of 0o to the total phase bias: 

d>1 = it + arcsin [ — ] , = n — arcsin ( — | . (12) 

\E J rx \E J 

(D FM (2) N FM 

In the FM configurations, the loop currents again cancel: (0 O ) = (0 O ) =0, while 
1 lVI are still given by Q. For the AF ones, 0^f are as in (jSJ) and (J7|) (though with a 
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different 4>\ i), yielding 



af ; 2(I p /2eE) 



(1/Eo) 



? 

Using these to expand the potential (fTUj) . one finds the proper generalization of 



(13) 



W AF ) = U 3 (d>™) - + O[{E/E f\ . (14) 

2eV^o - J 

Increasing the current bias decreases the effective Eq, mcreasing the coupling energy: as 
expected, a near-critically biased central junction is more sensitive to changes in current 
direction in the outer arms. Thus, tunability of the coupling strength (but only up) as in (j!4)l 
requires an extra current-bias lead for X itself, plus the generation of a flux-bias asymmetry 
as in ()12j) . While the latter is a common resource in two-flux-qubit experiments 0, E|, 
presently a substantial relative detuning of the coupling strength would require a flux-bias 
asymmetry of a significant fraction of a flux quantum. 
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In practical terms, the prediction is an interaction Hamiltonian H mt = Ja a z a h zl with 
coupling J = (h/2e)Ip/ a/ I^ —I^, in terms of the critical current of the shared junction I c0 
and the bias lb- Equivalently, the effective AF mutual inductance has absolute value M e g = 
(h/2e)/ a/^co — ^b) which again equals the Josephson inductance of the coupling junction at 
the working point. Linearizing these relations, we can sketch how an Jb-noise translates into 
a J-noise. One finds the width of the coupling energy as 

mill 

aj = nJ^^o) , 

4e 2 (/ c 2 -4 2 ) 3 

where the bias-noise power is S h (u) = Jdte iu3t {6I h (t)6I h (0)). This concludes the treatment 
of the coupling; further details depend on the circuit into which it is incorporated. 

Generalization to asymmetric devices and/or biases presents no trouble. It seems more 
interesting to further investigate coupling tunability. The standard ploy of a flux-biased 
compound central junction seems feasible, though it again does not directly lead to switchable 
coupling, which here would require E — > oo not E — > 0; also, flux leakage to the qubit 
loops is a possible problem. 

An attempt in this direction is presented in Fig. |21 The "qubit" junctions 1-6 are taken 
the same as before, while the "coupling" junctions 7-10 are all large, with Josephson en- 
ergy Eq. It is further assumed that no flux is threading the "crossbar" part of the circuit, 
that is, the legs which are drawn diagonally for clarity should actually lie almost on top 
of each other. This fixes 0io = 07 — 08 + 09 (fo r a convention in which positive phases 
correspond to a clockwise current in the a-loop). As before, also 02,5 are readily eliminated, 
while all stationary configurations have <fii = 03 and 04 = In terms of the remaining 
variables, the potential can be written as 

U(<f>) = -E[2 cos 0i + a cos(0"-20i-0 7 -0 9 ) + 2 cos 4 + a cos(0k-20 4 -0 8 +0 9 )] 

- E [cOS(f) 7 + COS 8 + COS 09 + COs(0 7 -0 8 +0 9 )] - Xi0 7 - J 4 8 + (J 3 +J 4 )0 9 . (16) 

The counterpart to (|T2*j) becomes that, to leading order in E/Eq, degenerate bias should 
be implemented as 




FIG. 2: A "crossbar" circuit in which 3JJ qubits are coupled with tunable strength and sign. 
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(15) 



The four bias leads enable a bewildering degree of tunability, subject to Ylj=i^-j = 0- 
However, only two one-parameter families will be considered, which are thought to be the 
most useful: (A) X\ = —I2 = —I3 — I4 — T, and (B) X\ = I2 = —I3 = — Z4 = I. In 
biasing scheme (A) one trivially finds 07°'* = = arcsin(X/£'o) and = 0: the external 
current flows through the horizontal legs. Conversely, in scheme (B), <fy^ = 0g ^ = and 
(j)^ = arcsin(X/£'o) so that the bias flows through the vertical legs. Minimizing U in (jl6j) . the 
leading corrections to these phases due to the qubit currents are readily found. Substituting 
these back into U, one obtains the energies of the various stationary configurations. 

Bias scheme (A) favours the FM states: 



U 



A,AF 



— u 



A,FM 



I 2 
P 



1 



4e 2 \^/El -J 2 



1 

Er 



(18) 



The interpretation is clear: the vertical, unbiased junctions mediate an AF interaction as in 
Fig. ^ The horizontal junctions mediate a "twisted" FM interaction, which overcomes the 
AF one due to the current bias. Compared to and (flljl. Eq. (JTSj) is a factor 2 smaller: 
the two junctions in parallel in Fig. El have a coupling energy which is twice as large, and 
therefore an effective inductance twice as small, as the single shared junction in Fig. ^ 

On the other hand, bias scheme (B) favours the AF states, with a coupling energy which 
is the exact opposite of the one in (|18j) . In fact, if one envisions flipping the right part of 
Fig. |21 it becomes clear that, when changing from scheme (A) to scheme (B), the role of FM 
and AF states is simply reversed. 




FIG. 3: Sign-tunable galvanic coupling of "black-box" qubits a and c via an intermediate fo-loop. 

Another scheme for sign-tunable coupling, more similar to the devices described in jT(| 
111 ], is shown in Fig. |3| We generalize and at the same time simplify the analysis by taking 
the two qubits as "black boxes", which may be asymmetric, biased away from degeneracy, 
and/or of a different design than 3 J J. This stresses that the crucial element in the analysis 
is the coupler, not the qubits. The potential reads 



U 3 (<j>) = C/ o (0x-0i) + U c (<f> c x +<f> 3 ) - E^osfa - £ 2 cos(#+0 a -0 3 ) - ^ 3 cos0 3 



(19) 



Here, the qubit potential U a is already minimized over its internal degrees of freedom; we 
denote U' a (<p£) = —I pa /2e and U"(<fi%) = Xa/4e 2 , where \a is the qubit susceptibility. These 
have two values depending on the qubit state, with I pa having two different signs; similarly 
for the c-qubit. The scheme is analogous to its magnetic counterpart in l3|, which means 
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-^-£ 2 sin<^ 


{ i ] = 




? ] = 


,(1) Xa / ,(1) 


s°= 





that i?2 is a variable of order one. We expand to second order in terms of only E 1 1 and 
E% , implementing the coupling per se. 

For this purpose, we write <p\ = (fyp + cj)^ 1 and 03 = </>g + </> 3 . Minimizing Z7j, one finds 

(20) 

^ J = ^ + ^sin^, (21) 

if^cos^, (22) 
^ = + (0i 1J -4 1} )^ cos# . (23) 

Substituting these into (jTUj) for Uj is straightforward. To C^-Efl) one finds trivial constants, 
terms oc I payC representing small shifts in effective bias, terms oc I pac representing small 
(Josephson) inductances, and some small corrections oc \a,c- Most interesting for us is the 
coupling 

^•■■ + tif'^"- < 24 > 

This is the expected result: expressing (9) in [l2| for the magnetic case as H [nt = 
MabMbcXfrlpalpc (which, incidentally, generalizes to the quantum case), (j23| follows un- 
der M ab l/Ae 2 Ei, M bc l/4e 2 £ 3 0, and Xb ^ 4e 2 £ 2 cos<^, the susceptibility of an 
rf-SQUID of vanishing inductance. Compared to the "crossbar" design, the proposal of 
Fig. |3] may be easier to implement. It also avoids realizing a small coupling as the difference 
of two larger ones, which may be sensitive to control errors. The restriction E% <C E\^ in 
the above can be lifted k|; the corresponding closed- form generalization of (|24|) will be 



presented elsewhere shortly. 

The derivations in this note have, hopefully, stressed transparent results and clarity rather 
than detailed modeling. For the latter, one should first of all study the classical potential 
nonperturbatively, since real coupling strengths will not be infinitesimal. One can {%) calcu- 
late the actual (coupling) energies at the codegeneracy point (symmetry point in parameter 
space) from a potential such as (JTJ), as is relevant to e.g. spectroscopy or at finite tempera- 
tures. Alternatively, one can {%%) calculate the classical stability diagram s(</» x ), indicating 
which flux state (s a , s&) = (1, — l]etc. prevails as a function of the biases; this is useful for, 
e.g., impedance measurement 0,111]. While (i) and (ii) are not equivalent outside the range 
where the energy-bias relation can be linearized, both are special cases of the classical "band 
structure" U(s, <fi x ), which is readily calculated whenever needed. 

Also, the classical derivation of (}2"4"|) etc. is a simplification. A full quantum analysis 
could be perturbative in the coupling or not, since its relevance depends on the junction 
capacitances, not on the classical interaction strength. Such an analysis can be transcribed 
from its magnetic counterpart ^2|, but this will not be pursued here. 

Extending the above designs to a linear qubit chain looks straightforward; a possible 
adaptation to 2D qubit lattices remains to be investigated. 
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